The probe readout and quantum limited measurements 
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Assuming that the parameter dependent evolution, as well as the measurements that are done for 
readout, of a quantum system that acts as the probe in a quantum limited measurement scheme are 
both fixed, we find the optimal initial states of the probe that will saturate the quantum Cramer- 
Rao bound. When the probe system is itself made of identical, elementary, two-level subsystems 
or qubits, we connect the optimal state of the N qubit probe to that of the one qubit probe. This 
is done for two different classes of dynamics for the probe qubits, one of which is entangling while 
the other is not. We study the limitations placed on the optimal initial state of the probe and the 
achievable measurement uncertainty by restrictions on the readout procedure that is applied on the 
probe qubits at the end of the measurement protocol. 



I. INTRODUCTION 



In modeling quantum limited measurements - in par- 
ticular single parameter estimation - a quantum system 
that acts as the probe interacts with the measured sys- 
tem in a way that depends on the value of the parameter 
of interest [lJ-Q ■ Knowing the change in the state of the 
probe that is generated by a parameter dependent Hamil- 
tonian, an estimate of the value of the parameter is ob- 
tained. However there is scope to fall into an infinite re- 
gression at this step in the modeling of the measurement 
because, after all, knowing the state of the probe means 
a measurement of the parameters that describe the state 
of the probe. So the question, what or who measures the 
probe, and how is immediately thrown up. Part of the 
answer lies in the observation that complete knowledge 
about the state of the quantum probe is not required 
to obtain the value of the measured parameter. Further- 
more the probe is typically assembled by putting together 
a large number of elementary quantum systems with low 
dimensional Hilbert spaces. For instance, without loss 
of generality, one can assume that the probe is made up 
of N two level quantum systems or qubits The re- 
sults of vonNeumann type projective measurements [f| 
on to the small, countable, set of basis states of the in- 
dividual probe units is sufficient to estimate the possibly 
irrational value of the parameter of interest to any degree 
of accuracy. The accuracy being dependent on the nature 
of the quantum probe, the number of elementary units 
in it and the number of times the probe is applied. To 
avoid confusion with the measurement of the parameter, 
we will refer to the measurements on the probe itself as 
a readout of the probe. 

The readout process is not fully understood either. 
This is evidenced by over a century worth of discussions 
and literature on topics ranging from the collapse of the 
wave function to decoherence and pointer states @, Q. 
However there is not much debate as to whether such 
readouts can be performed in the laboratory or not, since 
it indeed is routinely done. Even more sophisticated 
readouts that include POVMs can be done in the lab and 
they ultimately boil down to making a larger number of 
projective measurements on the probe 0, [l(| ■ 



The motivation for this Paper is the observation that 
not many types of readouts are possible given the avail- 
able technology in the design of a quantum limited 
metrology experiment. For instance, in virtually all pre- 
cision measurements using interferometers and light, the 
readout of the state of the light at the output ports is 
limited to photon counting [l], fllMTsI ] - There are sev- 
eral ways of using photon counting to accomplish read- 
outs of parameters, other than photon number, associ- 
ated with the state of the light in an interferometer like 
relative phase etc by combining known transformations 
of the state with photo counts. Homodyne and hetero- 
dyne measurements fl7l - [l9j on light are two examples of 
such techniques. In particle based metrology protocols 
like Ramsey interferometry also, the output state of the 
spins or atoms are subject to either Stern-Gerlach type 
or fluorescence type readouts [2(], [2l[ . 

Given the restrictions on the types of readouts that are 
possible, we find the optimal input state of a quantum 
probe made of many qubits for a fixed readout proce- 
dure and for different choices for the parameter depen- 
dent evolution of the probe. In the next section we review 
the quantum Cramer-Rao bound that forms the basis for 
the formulation of the problem. In section III the opti- 
mal input state for a non-entangling evolution is found. 
The optimal input state for an entangling evolution is 
discussed in Section IV followed by a discussion of our 
results in Section V. 



SATURATING THE QUANTUM 
CRAMER-RAO BOUND 



II. 



The quantum Cramer- Rao bound [22-25] gives the the- 
oretical lower bound on the measurement uncertainty in 
the estimate of a single parameter X as 



6X > 



1 



> 



1 



(i) 



where T is the quantum Fisher information. The sym- 
metric logarithmic derivative operator C is defined im- 
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plicitly by the equation, 



1 (r . r x dp x _ i 
-(£p x+px £) = — = Px . 



(2) 



The measurement uncertainty SX is quantified using the 
units corrected, root-mean-squared deviation of the esti- 
mate of X from its true value: 
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Let a general readout on the probe be described by 
a POVM with a one parameter family of elements -£?(£) 
such that 



J dCE(0 = 1- 



Let p(£,\X) = tr(E(£)p x ) be the measured probabilities 
for various outcomes of the POVM when the true value 
of the measured parameter is X. As shown in (2^, the 
quantum Fisher information is given by 



J- = max F, 
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where F is the classical Fisher information computed 
from the probability distribution for the measurement 
outcomes as 



F 



d£p(£\X) 



= / d£ 



dlnp(C\X) 
dX 

dp(£,\X) 
dX 



The maximization in the Eq. @ is over all possible read- 
out procedures (POVMs) on the probe. Such a max- 
imization is indeed a daunting task and even if it can 
be done, implementing the POVM that maximizes the 
Fisher information, thereby minimizing the measurement 
uncertainty, may, in all likelihood, be impossible to im- 
plement in the lab. The second inequality in ([l} cir- 
cumvents the maximization problem by placing an up- 
per bound on F in terms of the expectation value of the 
square of the symmetric logarithmic derivative operator 
£. This expectation value can be computed directly from 
the initial state of the probe and its parameter dependent 
dynamics, independent of the readout procedure. In the 
context of the current Paper it is worth reprising the se- 
quence of steps detailed in [22| that lead to this upper 
bound. We have 
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where we have used Eq. ([2]), the cyclic nature of the trace 
and the Hermiticity of L and E(£) to obtain the last 
equality in the equation above. Our focus is on the case 
where the second inequality in (fTJ is saturated. When 
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we have 
F = 
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tr 
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Using the Schwarz inequality, \ (x, y)\ 2 < (x,x) ■ (y,y) for 
the trace norm, we can write an upper bound on the 
classical Fisher information as, 



F < d£ti 
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t r (E(Z)C 2 p x ). (7) 
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The Schwarz inequality is saturated when 

for all £. Assuming Eq. (JS) holds, we get 

F I dt;tr(E(Z)£ 2 px) - tr(£ 2 px) = (C 2 ) 



(8) 
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and the second inequality in (jTJ) is saturated. Equa- 
tions (0 and © furnish the conditions on the the read- 
out procedure (POVM) such that a quantum probe in the 
initial state px that undergoes the parameter dependent 
evolution implicitly contained in £ will attain the quan- 
tum Cramer-Rao bound. Multiplying Eq. ((SJ by */ E(£) 
from the left and J px from the right we obtain, 
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a ) px = o, 



(10) 



for any p x . The above equation is satisfied if the readout 
is taken to be a set of orthogonal projectors, E(£), on to 
the complete set of orthonormal eigenstates of £. The 
are inverses of the eigenvalues of £ with 
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Condition ([S]) implies that are real. 

In situations where the readout procedure is fixed due 
to practical reasons or otherwise, one can now formulate 
the problem of finding the optimal initial state of the 
probe p x as follows. We limit ourselves to the rather 
common case where the readout is a complete set of or- 
thogonal projective measurements. Even if the readout is 
realized by a more general POVM, we assume that a suit- 
able Neumark extension [H, [l(| has been used to reduce 
it to complete set of orthonormal projectors. Given this 
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complete set of orthogonal projectors denoted as 

we can construct the symmetric logarithmic derivative 

operator corresponding to this readout procedure as 



The optimal initial state of the probe px for this readout 
satisfies the equation 
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(Cp x + Px C) = p' = -i[H, Px ], (12) 



assuming that the parameter dependent evolution of the 
probe is generated by the Hamiltonian, 



probe 



= XH. 



(13) 



In the following we will investigate the solutions of 
Eq. ([12")) for two choices of H, one that entangles the 
probe qubits and one that does not. The number of 
probe qubits N is the resource against which the per- 
formance of the measurement scheme is calibrated. The 
discussion here is quite general and can be applied to the 
case where the number of probe units itself is not the 
most important resource. For instance, in interferome- 
try with light the circulating power in the interferometer 
and not the number of photons is the crucial, limited 
resource [ll| . In other words, in our discussion, N is es- 
sentially a place holder for the relevant resource for each 
measurement scheme and a mapping between the real re- 
source and N can be found quite easily for most quantum 
limited metrology schemes. 



III. NON-ENTANGLING EVOLUTION OF THE 
PROBE QUBITS 

Let the parameter independent part of the Hamilto- 
nian in (TTBl that governs the time evolution of the ./V 
qubit probe have the form 



A' 



(14) 



where hy> is an operator that acts only on the j th qubit. 
By construction the time evolution generated by this 
Hamiltonian will not lead to entanglement between the 
probe qubits. Without loss of generality, using the free- 
dom to define a basis independently for the two dimen- 
sional Hilbert spaces of each of the individual qubits, we 
choose all the single qubit operators to be identical and 
equal to 



given above, we make the assumption that the read out 
procedure is limited due to practical considerations or 
otherwise to projective measurements along the states 
|+) and |— ) for each qubit where 
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(15) 



The symmetric logarithmic derivative operator for which 
this readout procedure saturates the quantum Cramer- 
Rao bound is then 
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where ji = { + ,—}. The optimal state of the N qubit 
probe corresponding to this dynamics and readout can 
now be found by solving Eq. (fl"2"|) . 

We first look at the case where the probe is made of a 
single qubit; i.e. N = 1. In this case 
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We write an arbitrary state of the probe as 

1/ 

Px = ^l 1 + a-iVi), 

where <7,, i — 1,2,3 are the Pauli matrices and the de- 
pendence of the state on X is hidden in the dependence 
of the coefficients a$ on the estimated parameter. Using 

|+)(+| = i(l+cr 1 ) and \-)(-\ = 1(1 - ai ), 

and the anti commutation relations of the Pauli matrices 
we have 



l sr x 1 

2 {£ > Px}= 4 



1 + ai 1 — <Zi 
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1 + ai 
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a7 + aT 
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a7 + ~+ 
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and 



i[H, px] 



(16) 
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Inserting Eqs. dT6j) and CEZl) into Eq. © wc find a solution 

as 
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a7 



= -a 2 , 



with 



ai = 0. 



(18) 



Once we choose to define the basis for the Hilbert space 
of each qubit so that the evolution Hamiltonian is as 



Since 03 does not appear in the right side of the equation, 
we are free to choose its value depending on the available 
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state preparation procedure for the probe qubit. Using 
Eq. (JT5J) we have 

£ = — Oal+X+I + aa|— ><— | = -a 2 a x . 

The quantum Fisher information is therefore 

I '= (£ 2 ) = tx{a%tp) = aj. 

The choice, a 2 = ±1 maximizes J- and positivity of the 
state of the probe qubit now requires that 03 = 0. So we 
find two possible optimal states of the one qubit probe 
for the given dynamics and readout as 
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Px = ± cr 2 ) 



(19) 



The results obtained so far are not anything new or 
unexpected. In Ramsey interferometers [2l|, [26| using 
atoms with two effective states in play forming qubits, 
the effective evolution of the probe units is modeled as 
rotations about the 03 axis in the Bloch sphere while 
the qubits themselves are initialized along the o\ or a 2 
directions. The optimal readout is then measurements 
on the individual probe qubits along a 2 or a\ directions 
respectively. The point of the preceding discussion is 
primarily to illustrate the means of obtaining the optimal 
state of the probe given that the readout is fixed. 

For dynamics generated by a non-entangling Hamil- 
tonian of the form given in Eq. (fT4")l , it is known that 
the Heisenberg limited scaling of 1/N is obtained for a 
"Schrodinger cat" state that is highly entangled 
However this assumes the ability to do a phase kick back 
operation after the parameter dependent evolution of the 
probe followed by projective measurements in order to 
implement, in effect, a readout on to a basis of entangled 
states. If we restrict the readout on each probe qubit to 
be along the basis given in (|15p , then we have to again use 
the approach discussed above to find the optimal initial 
state of a multi-qubit quantum probe. 

We briefly discuss the N — 2 case first. We have 
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and 



H (2) = ^{<Jz®t + t®a z ). 

Using the commutators and anti-commutators for tensor 

products of Pauli operators given in [27], and assuming 

(2) 

that the initial state of the two qubits, p x ' has the generic 
form 

( 2) I 

p x ' =-(101 + a t <Ti <g) 1 + bjl ® aj + CijOi <g> aj), 

we obtain sixteen algebraic equations (see Appendix [X] 
for details) from Eq. (TT21 by equating coefficients of 



corresponding operators. It is worth noting that the 
Schrodinger cat states, 
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(|00)±|11», 



with the corresponding density matrices 



(2) 1 

p± = - (1 <S> 1 + a\ <g> a\ ± cr 2 ® a 2 + a 3 <g> ct 3 ) , 

are not solutions of the equations we obtain. This again is 
symptomatic of the restriction on the readout procedure 
we have imposed. On the other hand it is straightforward 
to verify, as is done in Appendix [Al that 
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-(2) 
Px 



Px ® Px, 



is a solution, where px is the optimal state of the single 
qubit probe obtained in Eq ([T9")) . We find that 



£ = -2| 
so that 
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■2\ --)(-- I = -<Tl® 1-1®«71, 



T = (£ 2 ) = 2(1 ® 1 + (Ji <g> ax) = 2 = 4(A 2 i? (2) 
Generalizing to N qubits we again find that, 

P { P = P%\ 



is a solution of Eq. (fT2"| with 
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Detailed proofs of these results are given in Appen- 
dices |XT] and |X2j Significantly, we see that (£ 2 ) scales as 
N rather than as N 2 . So one does not reach the Heisen- 
berg limited scaling of 1 /N for the linear, non-entangling, 
parameter dependent dynamics of the quantum probe. 

The main point of the preceding discussion on a par- 
ticular example of non-entangling dynamics is to high- 
light the fact that with a restricted readout procedure, 
it might not be possible to go beyond the shot noise lim- 
ited scaling of 1/ yN for the measurement uncertainty 
even if the ability to initialize the quantum probe in ar- 
bitrary entangled quantum states is available. In other 
words, for implementing quantum limited measurements 
that beat the shot noise limit, we see that devising ways 
of doing possibly complicated readouts can be as impor- 
tant as control over the initial state of the quantum probe 
and its dynamics. 



IV. ENTANGLING DYNAMICS 

Now let us consider entangling dynamics for the probe 
qubits generated by 
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(20) 
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Note that this Hamiltonian does not belong to the family 
of non-linear Hamiltonians discussed in [28], [29| that leads 
to measurement schemes in which the uncertainty scales 
as 1/N k or l/7V fc ~ 1 / 2 with respect to N depending on 
whether the initial state of the probe is entangled or not. 

The readout procedure in this case is also the same as 
before with independent measurements of the individual 
qubits along the |±) axis. However, despite this restric- 
tion p x N is not, in general a solution to Eq. ([12")) . In fact 
for the entangling Hamiltonian in Eq. (|20l) one can show 
that (see Appendix [B]) if p x d ^ for d = 1,2,... is not a 
solution of the same equation giving the optimal state of 
a quantum probe made of 2c? (even number of) qubits. 

For N = 2, it is still worthwhile to find the optimal 
state even if it cannot be px <8> px ■ One possible solution 
for Eq. (|T2"j) is the state 

(2) 1 

P x = -(l®l + CllCTl<8)CTl+C23CT2®a'3 + C320-3®(72). (21) 

This state is pure when c\\ = C23 = c 32 = 1. How- 
ever for this state (£ 2 ) = 1, indicating that it does not 
saturate the quantum Cramer-Rao bound for the given 
dynamics and readout. In fact the restriction that the 
qubits are measured independently means that even with 
two qubits, the measurement uncertainty is not improved 
compared to the single qubit probe. The solution to 
Eq. (|12p obtained in (l2~Tj) is not unique either. For in- 
stance, another solution is obtained immediately from 
the state above by changing the signs of C23 and C32 which 

in turn swaps 1/A++ and 1/A For larger N one can 

solve the system of algebraic equations generated from 
Eq. (fT2j) by equating the coefficients of corresponding op- 
erators to find one or more optimal initial states of the 
probe that saturate the quantum Cramer Rao bound for 
the entangling dynamics. 

Rather than solving for the optimal state on a case- 
by-case basis, it is more useful to pursue solutions of the 
form p x N even though we have ruled out such solutions 
for even N. For odd N we can show that p x N is a solution 
of Eq. (fl2|) with the A's given by 



Ajij2...jjv AjjAj-j ■ ■ • \j N 

where 1/Xj, are the eigenvalues of C for the N = 1 case. 
The proof of this result is rather technical and long and is 
given in Appendix O This is not a very useful result be- 
cause 1/Aj, have values ±1 and so l/^iuV.-jjv = il for 
odd N. A simple computation shows that then C? = 1 
and so (C 2 ) = 1. In other words, no advantage is ob- 
tained in having N qubits in the probe rather than one 
if the N are initialized in the state p® . Even with the 
entangling evolution, in a real experiment, if the ability 
to initialize and readout the probe in states that are not 
simple tensor products is not available, then the reduc- 
tion in measurement uncertainty promised by quantum 
limited metrology is wiped out. 



V. CONCLUSION 

A general quantum limited measurement for estimat- 
ing a single parameter can be thought of having three 
stages. There is a preparation stage in which the quan- 
tum system that is acting as the probe of the measured 
parameter is initialized in a particular quantum state. 
The second stage is the parameter dependent evolution 
of the quantum probe and the last stage is the readout of 
the probe. The advantages of using specific - often entan- 
gled - initial states of the quantum probe was explored 
extensively during the initial phase of the development of 
the theory and implementation of quantum limited mea- 
surement schemes [J, [5|, I30M33I ] . How the dynamics influ- 
ences the measurement uncertainty was explored more 
recently gHH MM- 

This paper is focused on the third stage of a quantum 
metrology scheme when considerations, practical or oth- 
erwise, limit the types of readout that can be done on the 
quantum probe. This analysis is done in the limited con- 
text of a qubit based metrology schemes. Extensions to 
other quantum limited measurement schemes including 
interferometry with squeezed states, NOON states etc may 
also be considered. We see that arbitrary state prepara- 
tions and dynamics might not be particularly useful in 
delivering an improved measurement uncertainty if there 
are limitations on the readout. For instance, in the case 
of the N qubit probe evolving under a non-entangling 
Hamiltonian, the Schrodinger cat state turns out not to 
be the optimal state because of the restriction that the 
readout is limited to independent measurements on each 
of the N qubits. In the case of the entangling dynam- 
ics we see that when the initial state of the probe is a 
product state then with the same restriction as before on 
the readout, the performance of the measurement scheme 
is no better than what can be done with a single qubit 
probe. 
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Appendix A: Non-entangling Hamiltonian 

For two qubits, Eq. (fT2"j) leads to sixteen equations con- 
necting A++, A^ , A |_ and A__ to the fifteen coefficients 

<Zj, bj and Cy defining the state of the probe. These are 
obtained by equating the coefficients of operators of the 
form a a ® erg, a, ft = 0, . . . , 4 with do = 1- Using the 
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notation, 



K±±± 



1111 
± -: ± -: ± 



A ++ A + _ A_ + A — 

the equations we get are 

K +++ + K+—ax + IsT — | — 61 + K—+cxi = 
K-+- + K— + a x + K +++ bi + K+—Cu = -46 2 
+ ^+++ai + + = -4a 2 

if — + + if-+-ai + if+ — 61 + iT +++ cii = -4ci2 - 4c 2 i 
K +++ b 2 + K+— C12 = 4&i 

i^H 62 + if+++Ci2 = 4dl — 4C22 

K+++CL2 + X-+-C21 = 4ai 
isT-+_Oa + X+++C21 = 4cn - 4c 22 
if +++ a 3 + if_ + _c 3 i = 
if_ + _a 3 + K +++ c 3 i = -4c 32 
if+++&3 + K+ — C13 = 
JC+ — rb 3 + i\T +++ ci 3 = -4c 23 
A'+++c 22 - if — + c 33 = 4ci 2 + 4c 2 i 
K+++C23 + K — + c 32 = 4c i3 
K — +c 23 + K +++ c 32 = 4c 3 i 
K+++C33 - K__ + c 22 = 0. (Al) 

For the particular case in which 
(2) 

Px = Px® Px 

= i(l <g) 1 + CT 2 ® 1 + 1 ® er 2 + CT 2 ® CT 2 ), 

the non-trivial equations amongst Eqs. (|A1[) are 
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+++ 



K. 



0. 



= K+—= -4. 
From these equations we get 
1 1 



A-i 



A- 



-2, 



and 



1 



1 



0. 



A+- A — 1 
as a possible solution. Note that in this case, 

^ = ^ + ^, 3, k = {+,-}. 
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A, At 



1. N-qubits 

In [22j], an alternate expression for the symmetric log- 
arithmic derivative operator is obtained as 



£(o)= E 



Pj+Pk^O 



Pj + Pk 



O jk \j)(k\; O jk = (j\0\k), 



where O is the operator on which C acts. We write the 
optimal single qubit state as 

Px = ^2pj\j){j\- 

3 

In our particular example the basis in which px is diag- 
onal is given by the vectors, \i) — (|0) + i\l))/^/2 and 
V) = — l '|l))/v2- in the remainder we want to use 
the result J- = (C 2 (p')), also obtained in [22j. For the 
single qubit probe with evolution generated by H = er 3 /2 
we have 



Pjk = -HjWHi Px]\k) = --{p k -pj){j\a 3 \k), 



and 



-c(p') = -<E 



Pk - Pj 



k Pj +Pk 



0V3|*)|j)(*|, 



with the sum extending over all j, k such that Pj +p k =/= 
0. If the readout procedure corresponds to projective 
measurements with elements, 



we have 



tT[pE n C(p')] = 



Pk 



jk 



Pk ~ Pj 
Pk + Pj 



(j\a 3 \k)(k\0 n )(0 n \j) 



and 



tv[pE n ] = EftlO'l^ 



This gives us 

1_ £ jk Pk*5%ti\<r 3 \k)(k\0 n )(0 n \j) 

The N qubit state tensor product state can be written 

as 

Px = E p ^ p ^ ■ ■ ■ Pi" \ Jlj2 ■ ■ ■ JN ^ ( Jlj2 --'3n\, 

jl,j2—jN 

while the Hamiltonian corresponding to non-entangling 
evolution on the N qubits is, 
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®(N-1) 



From p' = — i[H, p], and denoting the string j\]2 ■•■3n 
as j, we have 

P y = -i{j\[H,p]\k) 



JV 



J^PnPn ■ ■ -Pk, ■ ■ -Pj N (jiW 3 \ki) Yl S . 



jnk n 



1 = 1 



N 



+ 9 E^' lP J' 2 • ■ - p i' ■ ■ - p i N Ossify) II ^" fe " 



JV 



/=1 
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using the above, we obtain 

N 

<w = -<E E 



Ph - Pji 
Pji + Ph 



(ji\°z\h) 



x \j 1 ,...ju...j N ){j u ...k u ...j N \. 

Assuming that the readout procedure consists of projec- 
tive measurements corresponding to the operators, 

-^n l^ni @ri2 ' ' ' @njv ) i^ni ' ' ' ®n N 1 3 



we have 



N 



tr[pE n £(p')] = -iJ2 E 



Pj 



Pji ~ Pk t 
Pji + Ph 



N 



= -*E E RPlJijmK 

i=l K {j m } m ^lm^l 



Pj 



Pji - Pk, 
Pji + Pk, 



E 

jiki 



where {j m }m^i stands for j\ ■ ■ ■ ,jl-i,jl+u ■ --Jn- Simi- 
larly we have 



Ui\O nt )(0 ni \ki)), 



tr{E H p)= E Y[p 3 J(jm\0„ 



31---3N 



x fE^i^i^) 



So the eigenvalues of for the N qubit probe under- 
going non-entangling evolution is 



1 N 



i=l 

iV 



St- 



(A2) 



/=i 



2. C for non-entangling dynamics 

Using the expression for l/A r -j from Eq. (|A2|) and the 
fact that corresponding to the optimal state of the one 
qubit probe, the eigenvalues of L are =Fl, we can write 
the symmetric logarithmic derivative operator on the TV 
qubit probe as 

C = -N\+... +)(+■■■ +\+...+N\ >( 1 (A3) 



Using |+)(+| = (1 + cri)/2 and |-)(-| = (1 - oi)/2 we 
get, 

1 - 

c = ~2N E^ - 2r ) [ 6{t + <?iT {N - r) ® (i - °iT r } , 

where r is the number of |— )(— | projectors in each term 
in Eq. (|A3|) and C is a shorthand indicating all terms that 
are tensor products of N — r factors of (1 + oi) and r 
factors of (1 — <ti). 

Once the sum is distributed over the tensor product, 
we get terms with N — q factors that are 1's and q ci's 
with q = 0, 1, . . . , N. We first focus on the sign of the 
various terms with fixed number of 1's and ci's in the 
expression for C. If out of the q factors of cti's, an odd 
number s of them come from the r (1 — 01) factors then 
the term, as a whole is negative. Now, out of the total q 
factors of <j\ we can pick s of them in q C s ways. Now, out 
of the N — q factors of 1 , r — s of them have to come from 
the (1 — txi) terms. These r — s factors can be picked in 
N ~ q C r - s ways. The remaining 1's and ai's come from 
the (1 + a{) terms. So the total number of ways in which 
one can construct a term with q <j\ 's starting from a term 
with r factors of (1 — <ti) is 

min(r,g) 

N/-y o \ s N — q/~i qr~i 



8=1,3,- 



We have for all q > 



JV-l 



JV 



^_JN - 2r) iv C r = -N2 N + 27V ^ 



r=0 

Hence 



N N min(r,q) 



£ =^EE E (N-2r)»-«C r - s «C s 

q=l r=0 s=l,3,"- 

For g = l the summation in the C reduces to 

N 

2^(/V-2r) lf ~ 1 C r - 1 = -2 N 

r=l 

Now for any q > 1, we have 

) 

(N - 2r) Ar -«CV- s 9 C s 

N min(g,r) 

= ^E E N ~ q Cr- S q C s 
r=l s=l,3, ■■■ 

N min(q,r) 

-2(/V-g)E E N - q - l e r - s -i q C s 

r=l s=l,3, ••• 



N min(q,r) 

E E 



N min(t/,r) 

-2 ? E E N ~ 9C rs q - l C s -i (A4) 



r=l 8=1,3, 
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To compute the sums in the equation above, we use the 
following results: For q > 1, £Li,3,- q C l = 2^ x and 

Ei=o,2,... q ~ 1(J i = 2Q ~ 2 - E q- <EU becomes 

AT min(q,r) 

£ £ (JV - 2r) "-«C r _. «C? a 

r=l s=l,3,--- 

^ E ,c - E ^ 

i=l,3,— 3=0, l,— 

q N-q 

-2(N-q) q °> E N ~ qC i 

i=l,3,— .7=0,1,— 
tf— 1 N-q-1 

-2q 9 ~ 1(J i E ^-^ = 

i=0,2, ••• j=0,l,-- 

Thus the symmetric logrithmic derivative for an N qubit 
probe evolving under a non entangling Hamiltonian is 



C = -{CI 1 



HN-l) 



real and so is the first part of each term in the double sum 
in ([Bl]). So the product term \[i{ji\(T3,\h)(h\0 ni )(0 ni 
has to be pure imaginary for 1/ to be real as required 
for saturating the bound on the Fisher information as 
discussed in Section. [TTJ However each term in this prod- 
uct comes form each qubit in the probe. So if we assume 
that each qubit in the probe is in the optimal state px 
corresponding to the N = 1 state, then for each qubit 
(j\cr 3 \k)(k\8 n )(d n \j) has to be pure imaginary so that 
again, the bound is saturated as assumed. This implies 
that when N is even then A^ are all pure imaginary and 
so the tensor product state p x N is not the optimal state 
of the probe corresponding to the entangling dynamics 
and readout procedure that we are considering when N 
is even. 



1. Optimal state of a two qubit probe 

Using the same notation as in Appendix [5] Eq. ([12 
reduces to the following sixteen algebraic equations, 



Appendix B: Entangling Hamiltonian 

When the parameter dependent evolution is generated 
by the entangling Hamiltonian, af N /2 and we consider 
an initial state of the probe of the form p x N , we have 



Now 



= -^YiPjiijiWalk) + lY[pj,(jiW3\k\ 



using the above, we obtain 



E 



UiPh +U1P31 



Iltokal*/) 



'x\ji,j2,---jN)(ki,k 2 ,...k N \, 



and 



nOikal*iX*il«»,)^n,lii>. (Bl) 



We also have 



(B2) 



tr(^p) = 2JJp Jl |(j I |e n| 



We get by dividing the right hand side of Eq. (|B1[) 
by that of Eq. (|B2|) . Note that the denominator in (|B2[) is 





+ K+ 


— ai - 


VK. 


K—+ 


+ K_ 


+ _ai - 


hK. 


K+— 


+ K+ 




YK. 


K-+- 


+ K- 


-+a\ - 


V K. 



K__ +Cll 



K+- 
K. 



+ 0-3 
-C22 



K 4 



-C\\ 
-C31 

-C33 
-C12 



_a 2 - 
K-+-a 3 - 
K +++ b 3 H 
K+— 63 - 
K+++c 33 - 

K++ + C21 
K++ + C12 

K +++ c 32 - 
K +++ c 23 H 



K-+-C21 

K+++C31 

K+ C13 

K+++C13 

K — +c 22 
- X_ + _a 2 
6 2 

K — +c 23 

^-- + C 3 2 






-4c 23 
-4c 32 



4c 3 i 
4ci 3 
-46 2 


-4a 2 




46i 

4ai. (B3) 



In this noted earlier it is easy to verify that 

p x = px <8> px does not give any solution for 1/A n . 
Solving these equations for the variables a,i,bj and Cij 
and applying in the general form of p for 2 qubits we get 
the optimal initial state for the two qubit probe. Multi- 
ple solutions are allowed and for maximizing (C 2 ) we look 
for pure state solutions. One pure state solution may be 



obtained by setting all a 

Cll 



bj and to zero except for 



1 and c 23 = c 32 = ±1 which then leads to 
1 1,11 



A, 



A__ 



= — 1 and 



A_^ 



= c, 



where c can be any real number, including 0. 
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Appendix C: 1/A for odd TV 

Starting from Eq. ((BiJ by that of Eq. flB2j| we find 
the eigenvalues of C(p') for an TV qubit probe undergoing 
entangling evolution as, 

■ v U. l P k i ( I\ l P"i-ll z Pn) 
1 _ l UnM} J! ,. 11 , 



^]\{h\^\ki){ki\6 ni ){e ni \^ 



(Cl) 



The optimal single qubit state px is diagonal in the eigen- 
basis \i)} of the o-i operator and so we can write it 

as, 

Px =Pi\i)(i\ +Pi\t)(i\, 

so that \ji), \ ki) = {\i), As before the readout oper- 
ators, \d ni )(d ni \ 's are given by |+)(+| and |— )(— |. Using 
the inner products, (i|±) = (l=pi)/2 , (i|±) = (l±i)/2, 
= —i/2, and }(— \i) = i/2 we get for any 
choice of 8 n 's, 



in} 1 



T^v (Pi +Pl) N 



(C2) 



Since (z | er 3 1 z) = { * | cr 3 1 z } = and (i\as\i) = 1, only 
terms with j = k contribute to the expression for pL-. 

jk 

and hence to L(p') and 1/A„ as well. Consider a single 
term in the double sum over {ji , fc; } in the numerator of 
equation (|C1|) . where 



and |fc) 



In this term we have the factor 



l[(3l\0n l )(en l \k l 



i N 



•T-±) 



In the double sum over {ji, k{\ in equation (|C1|) . there 
will be N Ci terms each having TV — 1 (=f) terms and one 
(±) term. In general, there will be N C r terms each hav- 
ing TV — r (=f) terms and r (±) terms. When Eft and 
hence 9 ni is fixed, then every term with TV — r =p terms 
and r ± terms has the same sign and so they can all 
be grouped together. We represent these terms that are 
grouped together as {=p • • ■ =p± • • • ±}, since the group is 



N—r 



labeled by r. Using this notation and Eq. (|C2[) . we have 



' n Pjv — Po N 



A ,7 



Pn+Po 



N 



C oP ° N {T T • • • T} 



+Q P \ r -\ P il 1 1 N C lP ^ 1 {TT---T±} 
Pn-i+Pi 



+ Q P ?- T , P lZ N CrPN-r{T ■ ■ ■ T ± • • • ±} 



PN-r + Pr 



,P/V-r Pr N/-~t N-rr, _i_ _,_ -,-1 



Pif-r+PP 



r, 1 

; — w^t c n-iP! l±±---±T| 

Pn-i+Pi 
^ q P^_P$_n CnP n {±± ... ±} 
Pn+Po 



where Q 



and p\ 



Pn- 



fer~ r (PlY- Now 



consider the first and last terms of the sum as a pair: 
Pn+Po 

-Q P f-4r N c NP Z{±±---±} (C3) 

P N + Po 

Since TV is odd, for a given choice of Eft, {=F T • • • T} and 
{± ± • • • ±} have opposite signs. Hence, equation (|C3[) 
reduces to 



A 1 = ±Q ( N C oP ° N - N C N p») 



Simillarly, we can pair up the remaining terms. Since for 
even and odd r 

Q ^- r ~/iZ N CrPN-A+- ■■ + + ■■■+} 
PN-r + Pr 

for any Eft's has opposite signs, we get 



1 _ ± i N+3 ( Pi -m) 



N 



Aft 



{Pi + Pi) 



N 



The sign of Aft depends on the choice of Efts. Using 
= +(Pi —Pl)/{Pi +Pl)i we see that for any odd TV, 



jN+3 



1 



Aft AnlAn2 ' ' ' A nN 



Generalizing the above result, the best initial state for 
any 2" (2m + 1) number of qubits can be obtained as 

p = pfn 2m , where m,n — 0,1,2,- •• and p<zn is the 
best initial state of 2™ qubits. 



10 



[1] V. Giovannetti, S. Lloyd, and L. Maccone, 

jScience 306, 1330 (2004)] 
[2] V. Giovannetti, S. Lloyd, and L. Maccone, Phys. Rev. 
Lett. 96, 010402 (2006). 

[3] C. M. Caves and A. Shaji, 

|Optics Comm. 283, 695 (2010) | 
[4] J. J. Bollinger, W. M. Itano, D. J. Wineland, and D. J. 

Heinzen, Phys. Rev. A 54, R4649 (1996). 
[5] S. F. Huelga, C. Macchiavello, T. Pellizzari, 
A. K. Ekert, M. B. Plenio and J. I. Cirac, 
|Phys. Rev. Lett. 79, 3865~(1 997) 
[6] J. v. Neumann, Mathematical Foundations of Quantum 

Mechanics (Princeton University Press, 1996). 
[7] J. A. Wheeler and W. H. Zurek, Quantum Theory and 

Measurement (Princeton University Press, 1983). 
[8] W. H. Zurek, |Rev~Mod. Phys 75, 715 (2003)] 
[9] M. A. Ne umark, Izv. Akad. Nauk SSSR 4, 277 (1940). 
[10] A. Peres, |Found. Phys. 20, 1441 (1990)| 
[11] B. P. Abbott et al, Rep. Prog. Phys. 72, 076901 (2009). 
[12] A. Chiruvelli and H. Lee, e-print arxiv:0901.4395 (2009). 
[13] U. Dorner, R. Demkowicz-Dobrzanski, B. J. Smith, J. S. 
Lundeen, W. Wasilewski, K. Banaszek, and I. A. Walm- 
sley, Phys. Rev. Lett. 102, 040403 (2009). 
[14] J. P. Dowling, Contemporary Physics 49, 125 (2008). 
[15] T.-W. Lee, S. D. Huver, H. Lee, L. Kaplan, S. B. Mc- 
Cracken, C. Min, D. B. Uskov, C. F. Wildfeuer, G. Vero- 
nis, and J. P. Dowling, arxiv:0909.3008 (2009). 
[16] B. C. Sanders and G. J. Milburn, 

|Phys. Rev. Lett. 75, 2944~(l 995) 
[17] P. Hariharan, Optical interferometry (Academic Press, 
1985). 

[18] M. Fox, Quantum optics : an introduction (Oxford Uni- 
versity Press, 2006). 



[19] K. J. Gasvik, Optical metrology, 3rd ed. (Wiley, 2002). 
[20] P. R. Berman, Atom interferometry (Academic Press,, 
1997). 

[21] S. Gleyzes, S. Kuhr, C. Guerlin, J. Bernu, S. Deleglise, 
U. Busk Hoff, M. Brune, J.-M. Raimond, and S. Haroche, 
|Nature 446, 297 (2007)] 
[22] S. L. Braunstem and C. M. Caves, 

|Phys. Rev. Lett. 72, 3439 (1994)| 
[23] S. L. Braunstein, C. M. C avei; and G. J. Milburn, 

|Ann. Phys. 247, 135 (1996)] 
[24] C. W. Helstrom, Quantum detection and estimation the- 
ory, 1st ed., Mathematics in science and engineering, Vol. 
123 (Academic Press, New York, 1976). 
[25] A. S. Holevo, Probabilistic and statistical aspects of quan- 
tum theory, 1st ed., North-Holland series in statistics and 
Probability theory, Vol. 1 (North-Holland, Amsterdam, 
1982). 

[26] S. Boixo, A. Datta, S. T. Flammia, A. Shaji, E. Bagan, 

and C. M. Caves, Phys. Re v. A 77, 012317 ( 2008). 
[27] C. Altafini, Phys. Rev. A 72, 012112 (2005)| 
[28] S. Boixo, S. T. Flammia, C. M. Caves, and J. Geremia, 

|Phys. Rev. Lett. 98, 090401 (2007)] 
[29] S. Boixo, A. Datt a, S. T. Flammia, A. Shaji, E. B agan, 

and C. M. Caves, |Phys. Rev. A 77, 012317 (2008)] 
[30] C. M. Caves, Phys. Rev. Lett. 45, 75 (1980). 
[31] C. M. Caves, Phys. Rev. D 23, 1693 (1981). 
[32] B. Yurke, Phys. Rev. Lett. 56, 1515 (1986). 
[33] B. Yurke, S. L. McCall, and J. R. Klauder, Phys. Rev. 

A 33, 4033 (1986). 
[34] S. Boixo, A. Datta, M. J. Davis, S. T. Flammia, A. Shaji, 

and C. M. Caves, Phys. Rev. Lett. 101, 040403 (2008). 



